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The electromagnetic mode spectrum of single layer graphene subjected to a quantizing magnetic 
field is computed taking into account intraband and interband contributions to the magneto-optical 
conductivity. We find that a sequence of weakly decaying gtiasi-transverse electric modes, separated 
by magnetoplasmon polariton modes, emerge due to the quantizing magnetic field. The character- 
istics of these modes are tuneable, by changing the magnetic field or the Fermi energy. 
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I. INTRODUCTION 

The field of plasmonics has been considered the pho- 
tonics milestone of the year 1998. To this choice con- 
tributed the landmark paper of Ebbesen et al. on the "ex- 
traordinary optical transmission through sub-wavelength 
hole-arrays". Jj| The effect was explained on the basis of 
the properties of surface-plasmon polaritons.Q 

Surface-plasmon polaritons (SPP) are electromagnetic 
surface- waves, guided by a metallic interface, resulting 
from the the coupling of the electromagnetic field to 
the collective plasma excitations of the metal. These 
guided modes are of important in fields as different as 
light guides at the nanoscale,[3[ spectroscopy and sensing, 
enhancement of light absorption in solar cells, enhanced 
Raman spectroscopy, and others. 

Graphene, being an one atom-thick metallic film is an 
obvious candidate for investigations on SPP. Recent re- 
search has established plasmon-bascd enhanced Raman 
spectroscopy and photocurrcnt,[4, 5] as well as room tem- 
perature prominent absorption peaks in the terahertz 
spectral range, [fil and nanoscopy of mid-infrared radia- 
tion confinement. [?[ 

Both in the two-dimensional (2D) electron gas and in 
graphene, the plasmon dispersion has a square-root de- 
pendence on the wave vector: ^2d{<1) oc The linear 
dispersion of the electrons in graphene, e{q) = ivpl, 
where vp is the Fermi velocity, implies that 512d(9) oc 
y/kpq, where kp oc ^Jn^ is the Fermi momentum and rig 
is the electronic density. Q If a grid of period L is super- 
imposed on graphene, plasmons of wave number q ^ 1/L 
can be excited. [2, 0] Furthermore, since f22D(Q) oc y/kpq, 
we expect the scaling relation ri2D(q) c>c nl^^ L^^^^ , 
which has been observed experimentally. @ 

Graphene has a number of advantages over other 
metallic thin films used in plasmonics; e.g., the ability 
of changing its carrier concentration using a gate, al- 
lowing a fine control over the frequency range for plas- 
monic excitations, and long propagation lengths as 
compared to conventional SPP. [lOj, Furthermore, in- 
homogeneous doping in a single graphene sheet allows 



the drawing of SPP propagation paths. Q 

When an external magnetic field perpendicular to a 
2D electron gas is applied, hybridization between cy- 
clotron excitations and plasmons occurs originating mag- 
netoplasmon modes. 12|,|l3| The presence of the magnetic 
field gives rise to strong absorption peaks making the dis- 
persion of electromagnetic modes very sensitive to the 
frequency. In the present work, we demonstrate that in 
addition to magnetoplasmon polaritons (MPP), graphene 
supports extremely weakly damped modes, which due 
to their resemblance to conventional transverse electric 
modes, are here referred to as gwasi-transverse electric 
(QTE). This article is organized as follows; in Sec. [ll] 
we overview the dispersion relation of electromagnetic 
modes supported by 2D electron systems. We revisit 
the simpler problem of zero external field, for which two 
types of modes can exist: SPP and weakly damped modes 
(transverse electric) with characteristics similar to pho- 
tons propagating in a dielectric. The magneto-optical 
response of graphene in the presence of a quantizing mag- 
netic field is described in Sec. IIIII The full mode disper- 
sion is calculated in the presence of disorder by employing 
the optical limit approximation to the conductivity. The 
losses, confinement, and polarization of the solutions are 
studied carefully. The outlook and conclusions are pre- 
sented in Sec. IIVI Technical details and derivations are 
given in appendices. 

II. DISPERSION RELATION 

We consider an infinite graphene film in the xy plane 
embedded in a dielectric medium of permittivity (per- 
meability) e (/i).[14] A static quantizing magnetic field is 
applied along the transverse (z) direction. We focus on 
electromagnetic modes propagating along the x axis. 
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The symbols have the usual meanings: w denotes the an- 
gular frequency, q is the complex longitudinal wavevector 
and K encodes the amount of confinement along the trans- 
verse direction. Maxwell equations relate these quanti- 
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ties according to the general relation = + eA*w^, 
so that in general both q and n are complex quantities. 
Note that a similar equation holds for the magnetic field 
H{r,t). Throughout, we employ SI units and the nota- 
tion z = z' + iz" for complex variables. 

The dispersion relation of electromagnetic modes fol- 
lows from the boundary conditions for the fields at the 
interface z = 0,[12] 
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where B is the intensity of the magnetic field, Z = -^//i/e 
is the impedance of the surrounding medium and ul (ch ) 
denotes the longitudinal (Hall) conductivity of graphene. 
The physical solutions of the above equation, uj = VL{q), 
contain the full mode spectrum of the system — a deriva- 
tion of the dispersion relation is given in Appendix A. 
The characteristics of the mode spectrum are determined 
by the conductivity tensor. The latter depends on 5, w, 
B, and, generally, also on the Fermi energy Ep, tem- 
perature and sample-specific broadening parameters. In 
the present work, we neglect the dependence of cfhh) 
on the in-plane wavevector g, and thus, hereafter, cr_f/(L) 
denotes the optical (local) limit of the dynamical conduc- 
tivity, i.e., (Th(l) = (^H(L){^,B) = aH{L)iO,^,B). The 
latter is justified for small wavevectors, more precisely 
for qls ^ 1, where Ib = y/h/eB denotes the magnetic 
length (— e < is the electron charge). 

In the absence of an external field B, the solutions of 
Eq. ([2|) are the so-called transverse electric and transverse 
magnetic modes, namely, k — zf2/itTL(J7, 0)/2 and n — 
2ifle/aL{^,0), respectively. In the former, the electric 
field is perpendicular to the direction of propagation, the 
mode dispersion is close to the light line, O ~ eg', and 
damping is small (g", k' <^ q'). Transverse electric modes 
require a negative reactive conductivity, a'[{u!,Q) < 0, 
and hence are not observed in conventional 2D gases (we 
note that when a'l < and (j'i^/\cr'l\ <^ 1 the behavior 
of the system resembles that of a dielectric). On the 
other hand, transverse magnetic waves are confined to the 
metallic surface, featuring large field localization, thus 
having important applications in sub-wavelength optics 
and plasmonics. (l5| In graphene both modes can exist 
because near at interband threshold, uj ~ 2Ep/h^ the 
function a'l{uj,Q) changes sign.[l^ Coupling and guiding 
of transverse electric modes by graphene in zero- field have 
been recently reported in Ref. |17 . 



III. MAGNETO-OPTICAL MODES 

When a magnetic field is turned on, electrons acquire 
considerable cyclotronic energies via the Lorentz force, 
and at sufficiently high fields, the continuum Dirac quasi- 
particle spectrum condensates into Landau levels (LLs) 
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Figure 1: Doped graphene in a magnetic field field of 5 T at 
zero temperature [Ef = 0.1 eV). Left: ctl is plotted as a 
function of the frequency / = u}/2ti. Clearly, a'l <Q in vari- 
ous frequency intervals. The full quantum calculation is seen 
to be crucial above ~ 20 THz. Right: The first few electronic 
transitions are shown. The Fermi energy lies slightly above 
the Landau level with n = 1, and thus only transitions to lev- 
els with n > 2 are allowed. The (small) dashed arrow stands 
for the intraband transition responsible for the strong peak 
observed near at 7.8 THz. Other arrows represent interband 
transitions. 



[Fig. [T] (right panel)]. In order to investigate how the 
zero-field picture changes due to the quantizing mag- 
netic field, we employ linear-response theory within the 
Dirac cone approximation to obtain an expression for 
'^L{H)(j-^T B) with both intraband and interband contri- 
butions included. (See Appendix B.) 

The main features of aL{uJ,B) for doped graphene in 
a quantizing field are shown in Fig. [1] (left panel) for fre- 
quencies up to the near-infrared. As a figure of merit, 
we choose Ep = 0.1 eV and field intensity i? = 5 T, 
both being typical values used in experiments. aL{uj, B) 
consists of an intraband term with spectral weight lo- 
cated at the lower-end of the spectrum (the strong peak 
at / ~ 8 THz), and interband high-frequency contribu- 
tions originating a series of peaks above the interband 
threshold / ~ 43 THz. The position and shape of each 
of these peaks depend on the particular values of Ep and 
B, but the main qualitative features in Fig. [T] are general. 
For a comprehensive discussion of the magneto-optical 
response of graphene see Refs. 3, 1^. 

In order to account for disorder, we have used an 
energy-independent LL broadening of F = 2.5 meV (and 
also F = 10 meV in one figure) in consistency with the 
values found in pump-probe experiments performed in 
epitaxial and exfoliated graphene samples. pc" 
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on infrared spectroscopy studies of the Drude conductiv- 
ity of graphene. 22| The renormalization of the optical 
conductivity due to the electron-phonon interactions is 
neglected. From the studies of magneto density of states 
taking into account the optical mode at 200 meV.|23| 
we expect the latter approximation to be valid in the en- 
tire range of frequencies considered here. 

In 2D electron gases, plasmons and cyclotronic excita- 
tions (with frequency ujc) hybridize leading to the well- 
known semi-classical magnetoplasmon spectrum, 51^ = 
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12l . |l3| In view of the strong contribution of 



100 



interband transitions to a'l (e.g., see the discrepancy be- 
tween the semi-classical calculation and the quantum 
formula, even at low frequencies, in the left panel of 
Fig. [T]), this formula should should be of limited appli- 
cability in graphene. Moreover, Shubnikov-de Haas os- 
cillations originate many frequency regions with a'[ < Q, 
and hence, similarly to 2D magnetized electron Fermi 
we may expect the splitting of the mode spec- 
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trum into many branches. 

In what follows, unless stated otherwise, we refer to 
the particular system in Fig. [T] We complement the 
discussion with analytic expressions that can be used 
to compute the several quantities for general parame- 
ters. The full mode spectrum is presented in Fig. [2] only 
non-exponentially growing solutions of Eq. ([2|) are shown 
(k' > 0). The most notorious feature is the existence 
of a series of well-defined branches, labeled by the inte- 
gers n — 1,2, etc. These branches can be divided into 
two distinct sets, namely, the set of branches with disper- 
sion close to the light line (n odd) and the remaining (n 
even) . The former will be shown to have the basic prop- 
erties of transverse electric modes and hence are termed 
QTE, whereas the latter are MPP modes (with polariza- 
tion not necessarily similar to transverse magnetic modes 
due to the applied magnetic field). 

Our results borne out two peculiar features of 
graphene: i) n even branches have two distinct solutions 
for each wavevector q'. This degeneracy is a result of 
hybridization between even and odd modes (see below), 
and ii) the frequency domain size of each branch is non- 
uniform due to the structure of LLs in graphene. The 
first branch occupies a region [0, the second [/i,/2[, 
etc., where /„ is defined to be the n-th node of the re- 
active longitudinal conductivity, (T^(27r/„, B) = 0. For 
the system under discussion, the nodes read /i — 8 THz, 
/2 ~ 39 THz, etc. 

Let us know discuss with detail the region / < 20 THz 
spanning two branches, n = 1,2. Here, the magneto- 
optical transport is predominantly semi-classical (see 
Fig. [1]), and hence the relevant frequency scale is the 
cyclotron frequency, uJc = evpB/\Ep\ {vp ~ 10^ m/s 
is Fermi velocity of carriers in graphene). At / < /i, 
the dispersion curve is pinned to the light dispersion line, 
f2(g') ~ cq' , except for frequencies very close to /i, where 
the dispersion curve shifts slightly towards the left-hand 
side (see inset) and part of the transverse electric char- 
acter is lost. The latter signals the onset of a rapid in- 
crease of a'l, as a result of a strong absorption peak at 
/i ~ a;c/27r (Fig. [1]). Last, when / becomes larger than 
fi, a'l changes sign, and a MPP branch develops. This 
mode displays the strong confinement to the 2D inter- 
face {k' ~ q') that characterizes conventional transverse 
magnetic modes. Its polarization varies significantly with 
frequency /wavevector as shown later. 

The simultaneous presence of the two distinct branches 
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Figure 2: Mode frequency / — Q/2n for doped graphene in 
a field of 5 T is plotted as a function of the longitudinal mo- 
mentum q' [solid (black) line]. For completeness, the light 
dispersion, / = cq' /2n is shown by the dashed (red) line, and 
the semi-classical solution is given by the dashed-dotted (blue) 
line. Inset: Mode spectrum near zero momentum. Arrows in 
the main panel indicate QTE solutions and numbers identify 
the distinct branches. Other parameters as in Fig. \T\ 



{n = 1, 2) in the intraband region can only occur for suffi- 
ciently high fields, B>Bc = \Ep\r/ehvp, otherwise one 
obtains a a single type of solution with plasmon char- 
acter, fl{q') ^ ^/q', as for B = 0. 25, 2^ The reason is 
that for B < Be the reactive part of the conductivity 



is always positive below the interband threshold, forbid- 
ding the existence of modes with dispersion close to the 
light line (QTE modes) in the low-frequency region. The 
transition frequency (at which the MPP n — 2 mode 
begins) is given by the solution of Eq. ^ with q' = O.pTj 



(3) 



Using the above expression, we obtain f2(0) w 27r • 8 THz 
in good agreement with the numerical result (see inset in 
Fig.©. 

The semi-classical expression for the MPP dispersion 
f2(MPP) (9') can be derived assuming T — and ignor- 
ing the interband contributions to the magneto-optical 
conductivity. 
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where il2B{q') = (e/?i) q'\Ep \ / (27re) is the graphene's 
longitudinal plasmon dispersion in zero field. The above 
expression is valid for q' 3> ^/c, and for $7 within the 
regime of validity of the semi-classical transport theory, 
il <C 2Ep/h (see dotted-dashed line in the main panel of 
Fig. (2). Eq. predicts an increase of MPP's frequency 
due to the presence of a magnetic field. A derivation of 
the above formula is given in Appendix C. This result 
coincides semi-classical magnetoplasmon spectrum for a 
2D electron gas.[i3,|ll| 

We now turn our attention to the mid/near- infrared 
part of the spectrum, where new branches (n = 3,4, ...) 



4 



emerge due to interband transitions. Figure [5] shows 
that the quantum calculation (solid line) deviates con- 
siderably from the semi-classical result already at / w 
20 THz. In particular, the quantum corrections cause 
a considerable slow down of the MPP's group velocity, 
Vg = dq'^l{q'), relative to its semi-classical value. This 
effect comes from the superposition of interband reso- 
nances tails that contribute with substantial weight even 
well below the interband threshold. For instance, in the 
range / « [20, 30] THz, the interband terms yield a cor- 
rection to the conductivity of about ~ 0.5i e'^/h (see 
Fig. [T]) explaining the bending of the solid curve rela- 
tively to the dashed-dotted curve in Fig. [21 Near at 
/ = /2 ~ 40 THz, a'l changes sign again, and a large 
bandwidth (wlO THz) QTE mode develops. The first 
six QTE modes are indicated by arrows in Fig. [2] Their 
dispersion relation is well approximated by ^{q') — cq' , 
except within the QTE/MPP crossovers (/ ~ /„, with n 
odd), where ^{q') acquires a complex form due the strong 
variation of the optical properties induced by sharp ab- 
sorption peaks in these regions. 

The full dispersion relation for the MPP branches is 
rather cumbersome because, as noted above, many inter- 
band terms contribute to the spectral weight around a 
particular frequency; see Eq. (|20p and text therein. A 
compact expression for g(w) valid for every MPP branch 
can still be obtained by considering T = and neglecting 
the Hall conductivity term in Eq. ^ . These approxima- 
tions are justified since i) for a quantizing magnetic field, 
the conductivity of graphene does not vary significantly 
with temperature, and ii) aH can be shown to provide 
a small correction only in the vicinity of each /„. We 
obtain. 
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In these expressions, E„{B) — Sn^/^\n\hvF/lB [where 
Sn = sign(n)] is the energy of the n-th (degenerate) LL 
in graphene, Np = mt[{EF / Ei{B)y] > is the number 
of occupied(empty) electron(hole) degenerate LLs with 
\En\ > 0, and A„(i?) is the n-th interband resonance en- 
ergy, defined as A„(B) = En+i{B) + En{B). The prime 
in the summation sign indicates that if Np ^ the first 
term is to be halved. Also, a cutoff n < Ncut must be 
taken when computing this summation (see Appendix 
B). For simplicity, the above expression for 4'(i3,w) only 
includes the interband contribution to ctl . The inclusion 
of the intraband spectral weight [see Eq. ([M)) ] in '^'{B, lj) 
is straightforward and plays a role only in the first MPP 
branch. Figure. [3] shows that the MPP spectrum com- 
puted from Eq. ([5]) cannot be distinguished from the full 
calculation. The latter agreement extends down to low 
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Figure 3: Mode frequency for doped graphene {Ef = 0.05 eV) 
in a field of 6 T as a function of q' [solid (black) line]. The 
broadening reads F = 10 meV. Arrows in the left panel indi- 
cate the first three QTE modes. 



frequency (the first MPP branch) because in the system 
plotted in Fig. [3] all electrons (and holes) occupy the zero 
energy state, that is, Np = 0. Indeed, according to our 
definition of interband and intraband contributions to the 
conductivity when Np = (see Appendix B), '^{B,uj) al- 
ready contains the full spectral weight in this case. 

We have discussed the general features of the mode 
spectrum of graphene under a quantizing electromagnetic 
field. It has been shown to consist of several branches, 
each possessing two distinct types of modes. In what 
follows, we demonstrate that MPP solutions have con- 
ventional decaying properties of SPP, whereas the QTE 
modes are essentially non-decaying, with electric field 
nearly transverse (hence their name). The MPP solu- 
tions will be shown to have a rich polarization diagram 
with longitudinal (transverse magnetic) character only 
emerging at specific frequency intervals. 



Decaying and polarization properties 

The decaying properties of the modes are summarized 
in Fig. [31 QTE modes display large localization length in 
X [z] direction, namely, q" /q' [k' /q'] in the range 10^^- 
10~^ [10~^-10~^] (sec for instance the first branch; top 
left panel). MPP modes, on the other hand, always show 
considerable decay along the z direction, k' = 0{q'), indi- 
cating strong confinement. The losses in the propagation 
direction x, on the other hand, vary appreciably and are 
determined by the graphene's absorption at the specific 
MPP frequency; e.g., in the range 8-20 THz, the longitu- 
dinal decay rate q" varies in the range O.lg'-lOOg', with 
maximum loss occurring near at the cyclotron frequency 
(~8 THz). 

Although the losses and confinement reported here 
have orders of magnitude comparable to those in the ab- 
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Figure 4: The ratios g"/?' and k' /q' , which characterize the 
longitudinal and transverse decaying rates, respectively, are 
plotted in logarithmic scale as a function of the mode fre- 
quency. Several panels are given to help visualize the differ- 
ence between QTE and MPP modes. (System parameters as 
in Fig.m) 

sence of a magnetic field [lo|, the strong dependence of 
these quantities on the frequency is exclusive to the 2D 
interface subjected to a strong external magnetic field 
(Fig. m shows that the decay characteristics can vary by 
several orders of magnitude around at / = /„ for all 
n). An important effect of the magnetic field is to al- 
low for QTE modes with lower losses than the zero-field 
transverse electric mode in specific frequency intervals; 
for instance, above 50 THz, q" / q' can reach a minimum 
value of the order of 10~^ , whereas for _B = (and other 
parameters unchanged) this ratio is about ^ 10^^ in the 
frequency window 50-100 THz. 

In order to complete our study, we demonstrate that 
the electric field of QTE modes are essentially transverse 
and study how the longitudinal (transverse magnetic) 
character of MPP modes depend on the wavevector. To 
this end, we compute the ratios £xy(zy) = Elx(z)/E!y and 
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Bx(z)/By. From Maxwell equations, we easily 
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and £{B)zy — £{B)xy x iq/K, where it is assumed B > 
so that an ^ 0. These quantities are plotted in the 
bottom panel of Fig. |4] for a frequency range spanning 
the n = 1 and n — 2 branches. Below 8 THz, the electric 
components ratios £xy(zy) (plotted in the right panel) are 
found to have magnitude in the range lO^'^-lO"^ (10~^- 
10"^), confirming that the electric field of QTE modes 
lies prominently along the y axis, resembling pure trans- 
verse electric modes [which have £xy{zy) = 0]. The mag- 
netic ratios of these modes are shown in the left panel. 
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Figure 5: Polarization properties of the electromagnetic 
modes within the first two branches. The left panel shows 
the complex modulus of the magnetic ratios Bxy and Bzy, 
whereas the respective electric counterpart, £xy and £zy, are 
given in the right panel. (System parameters as in Fig. [T]) 



As for the MPP branch, we find magnetic ratios 
Bzy in the range 0.5-1, and hence the polariza- 
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tion of these modes are distinct from transverse mag- 
netic modes (which have Bxy(zy) =0). In this MPP 
branch, the magnetic ratio Bxy decreases with increas- 
ing frequency/g' with a minimum (not shown) occurring 
at 34 THz, Bxy — 0.05. The latter signals the beginning 
of a crossover to a QTE branch: the longitudinal charac- 
ter is lost (Bxy increases) as the turning point (maximum 
q') is approached. Similar conclusions can be drawn for 
the remaining branches. 

The situation described in the above paragraph is dis- 
tinct from 2D electron gases, 13, 24 1 where solutions be- 
come more longitudinal with increasing wavevector. A 
question that deserves further investigation is whether 
the effect of finite q = q'sx in crj;, can influence the solu- 
tions at large q'. Note that Ib = O{10^^B^^/^), with 
B in Tesla, and hence a considerable renormalization 
of spectrum is expected for the largest wavevectors dis- 
cussed here, specially in Fig. [2] where q' can be as large as 
8 X 10^ m"^ The lack of longitudinal character of MPP 
at large q' reported here may indeed result from an inade- 
quacy of the optical limit in describing that region of the 
spectrum. On the other hand, the system represented 
in Fig. El (B = 6 T and r = 10 meV) displays modes 
with smaller wavevectors, making the optical limit less 
restrictive in this case. Notwithstanding, similar features 
are observed in this system, thus providing further evi- 
dence for the generality of the phenomenona discussed in 
this article (QTE-MPP crossovers, degeneracy of MPP 
branches, etc.). 



MPP wave localization 

We briefly address the wave localization characteris- 
tics of the MPP waves reported here. It is a well es- 
tablished fact that SPP in a metal can have wavelengths 
considerably smaller than electromagnetic waves of the 
same frequency in a dielectric. [l5j| In graphene in zero 



6 



1000 



500 




B=0 



20 



40 60 
f (THz) 



80 



11 



100 



Figure 6: The wave localization ratio A is plotted as a func- 
tion of the mode's frequency for a magnetic field of 5 T. The 
dashed-dotted (red) line shows the result for 5 = for com- 
parison. Other parameters as in Fig. [T] 



field this shrinkage effect is enhanced when compared to 
conventional 2D electron gas SPP.dSim Figure El shows 
the ratio of the wavelength in vacuum to the MPP mode 
wavelength, A = Aq/A (here, Aq = 27rc/w and A — 27r/q'). 
It can be seen that for 5 T, below the interband thresh- 
old, A is comparable to that obtained in the case of zero 
field. In particular, near the frequency resonant to the 
n = —1 n = 2 interband transition (see also Fig. [IJ, 
we obtain a large peak of about A ~ 10^, a figure compa- 
rable to the B — scenario (see how the peaks compare 
in Fig.ini). The remaining MPP modes in that figure show 
peaks with A ^ 10^. QTE solutions, on the other hand, 
have A ~ 1 regardless their frequency, a characteristic of 
transverse electric modes. 

A simple formula for A can be obtained in the intra- 
band region; using Eq. Q we obtain 
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where ag = e^/(47refi,c) denotes the graphene's fine struc- 
ture constant and A(i3) = hevpB/Ep is the cyclotron 
energy. The above expression predicts a decrease of A 
with the magnetic field, which is consistent with the 
curves reported in Fig. |5] in the region where Eq. ^ 
[and hence Eq. (O] is valid, i.e., 8-20 THz. In the limit of 
B ^ 0, Eq. ([9]) reproduces the result reported in Ref. 11 1. 



IV. OUTLOOK AND CONCLUSIONS 

We have computed the mode spectrum of electromag- 
netic modes supported by a graphene interface in the 
presence of a quantizing magnetic field. We have found 
a rich structure with extended crossovers between quasi- 
transverse electric (QTE) and magnetoplasmon polari- 
ton (MPP) modes as a consequence of characteristic 
Shubnikov-de Haas oscillations in the magneto-optical 
response of graphene. Analogously to the 2D Fermi 
gas, the dispersion relation splits in a large number of 



branches. 12|, |24| Interband transitions between the bot- 
tom and top Dirac cones originate terms in the conductiv- 
ity with considerable spectral weight in the semi-classical 
(low-frequency) region. As a consequence, the conven- 
tional semi-classical 2D magnetoplasmon dispersion be- 
comes restricted to a narrow wavevector / frequency inter- 
val. Our calculation within the optical limit approxima- 
tion to the conductivity predicts that, unlike 2D Fermi 
gases, each MPP mode with a given wavevector admits 
two possible values of frequency. The consequences of a 
non-zero wavevector q in the conductivity for the QTE- 
MPP spectrum is a challenging question and deserves 
further investigation. 

In summary, we have shown that a quantizing magnetic 
field changes the conventional picture of electromagnetic 
modes in graphene. 13, iBj In the LL regime, we predicted 
a series of branches in the mode spectrum, consisting of 
magnetoplasmon polaritons and quasi transverse electric 
modes. At small frequencies, a QTE persists, even in the 
semi-classical regime, as long as the cyclotronic frequency 
is larger than the electrons relaxation's rate. Due to the 
high LL energy gaps, these effects should be observable 
up to room temperature. 
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VI. APPENDIX A: MODE SPECTRUM 

For our purposes it is sufficient to consider a single 
Fourier component of the electromagnetic field in the fol- 
lowing form 



(10) 



where Km accounts for possible attenuation in the trans- 
verse direction to graphene. The subscript m = 1(2) 
denotes the region of space with z > 0(z < 0). Here, we 
consider the more general case of graphene embedded in 
dielectric mediums with permittivities ei and £2- 

We require > and q'^/Qm — 0- The former means 
that the wave can be confined to the z plane and must not 
diverge as z — >■ ±00, whereas the second condition allows 
wave attenuation as it propagates in the graphene plane 
(the formal divergence at x — — 00 is a consequence of 
the beginning of the perturbation at t = ~oo). 

Macroscopic Maxwell equations imply the following re- 
lation between the field amplitudes. 



(11) 
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where J-^ is the six-dimensional vector of amphtudes 
7'm = {Em,Bjn) and, 



Mrn,R 
M^r 



with 



1 



M, 



m.R 



(-1)'"^™ 

iqr,, 

-ig„i 

(-1)™+!^™ 

~iq„ 
ig,„ 



(12) 

, (13) 
(14) 



A straightforward consequence of Eq. pT|) is the weh- 
known relation, 



The mode spectrum for this problem is derived by impos- 
ing the boundary conditions for the electromagnetic field 
at the interface z — and making use of the relations 
between the field components [Eqs. (fTT|) - dTil) ]. The con- 
tinuity of the tangential (normal) component of the elec- 
tric field (magnetic induction) imply that qi = q2 = q, 
I^i,x{y) = ^2,a;(i/) ^ud Bi ^ = i?2.z- The discontinuity of 
the tangential component of the magnetic field yields 

— (^XxEl.x — (TxyEl,y , (16) 



1 1 

Bl.y B2^y 

Ml M2 



1 1 _ 

Bi.x B2,x — O-yxEi^x + '^yyBl.y , (17) 

Ml M2 

where CTq^ = aai3{q,uj) [a,/3 = x,y] denotes the dynami- 
cal conductivity of graphene. 



CTQ^(q,w) 



(18) 



and relates the Fourier transforms of the surface current, 
Js,a{qi^), and that of the electric field, Es^p{q,uj), at 
z — 0. The use of the local limit of the conductivity 
'^L{H){q = 0,aj) is justified whenever the wave vectors 
of interest are much smaller than the inverse of typical 
length scales. In the absence of a magnetic field, the local 
limit is justified for \q\ <^ , with aq being the carbon- 
carbon distance in graphene, whereas for a quantizing 
magnetic field the condition changes to |q| <C where 
Ib denotes the magnetic length (see Appendix B). 

Combining the above results, it is straightforward to 
obtain the general dispersion relation 



iuj 


[(- 


Hi £2' 








— ei 




)- 


f^L- 


X 






K2 ei, 




lUJ€l 














1 " 
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h-( 














^Mi 


M2 / if^J 



'H 



(19) 



In order to obtain the above expression, we have invoked 
rotational symmetry and use the notation employed in 
the main text: = cTxx and uh = <Jxy Setting 
ei(Mi) = £2(^2) = ^(m) (and hence ki = K2 = k) leads to 
the Eq. ^ in the main text. We remark that the term 
with an is negligible for most choices of parameters. We 
have verified that only near at the QTE/MPP transi- 
tions, where a'l ~ 0, the Hall conductivity an provides 
a small correction to the spectrum. 



VII. APPENDIX B: MAGNETO-OPTICAL 
CONDUCTIVITY OF GRAPHENE 

Within the Dirac-cone approximation, Q and model- 
ing the effect of disorder by an energy broadening func- 
tion, the magneto-optical conductivity of graphene at 
Fermi energy Ep and temperature T assumes the sim- 



(15) 

pie form[l 



<^L{H){^^,B) ^ gs9v ^ 

^ -lTh) npiEn) - upjE^) ^^q) 

where gs(v) = 2 is the spin (valley) degeneracy factor of 
graphene, npiE) 1/[1 + e(^--^^)/'=s^] stands for the 
Fermi distribution function, r„,„(w) is the LL broaden- 
ing, Anm = En - Em, with LL energies En given by 



E„ 



sign(n) [hvp/h 



(21) 



with Ib denoting the magnetic length, Ib = y''h/(eB), 
vp ~ lO^m/s is the Fermi velocity, and 



-(1 + 5mfi + 5nfi)5\ m\ — \n\ ,±1 ; 

(22) 

iS2'"((5|,„|J„|_i - <5|.m|-l,|ri|) ■ (23) 



/2 



The use of the low-energy (Dirac-cone approximation) 
theory to compute <yL{H) assumes an infinite sea of neg- 
ative energy states, and thus requires a cutoff |m| < 
iVcut in Eq. (PU)) . The respective cutoff energy E^^^t is of 
the order of graphene's bandwidth. Results are largely in- 
sensitive to the precise value chosen for the cutoff; in our 
numerical calculations we have considered En^^^ = 3 eV. 

The magneto-optical conductivity of graphene has two 
types of terms: i) intraband contributions corresponding 
to transitions within the same Dirac cone (i.e., n = m ± 
1), and ii) interband transitions that couple the valence 
and conduction Dirac cones (i.e., n — ~m ± 1). Transi- 
tions involving the zero energy LL (e.g., n = — > n = 1) 
need to be considered separately because this LL state 
contains electrons and holes. Here, for convenience, we 
classify the transitions involving the zero energy LL as 
interband-like. 
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The general expression Eq. (PO)) can be put into a more 
useful form by separating interband and intraband con- 
tributions (with the proviso made in the previous para- 
graph). For sake of simplicity, we assume T = and 
Ep > [the conductivity for holes can be obtained 
using the symmetry relations: (Jl{—Ep) — UL^Ep) 
and (Jh(—Ep) = —(Th{Ep)\. We denote the number 
of occupied electron LLs with > by Np, that 
is, Np = int \{Ef/Ei)'^'\ > 0. Intraband transitions 
(n = Np — > 71 = Np + 1) involve an energy difference 
of Aintra = ^/2hv p / 1 b{\/ N p -\- 1 — y/Np). Its Contribu- 
tion to the conductivity dominates at small frequencies 
where most of the spectral weight is concentrated around 
uj = Aintra/^- The intraband (semi-classical) conductiv- 
ity consists of a single term in the summation Eq. ((20)) 
and reads 



(intra) 



(intra) 



2e' 



1 — ihuj/T 



7^ r A 

^mtra 



2 fc2„,2 



[i-ihuj/ry 



ALtra/r^ 



ilr^ (1 



ihuj/Ty 



ALra/r^ 



(24) 
(25) 



According to our classification the latter equations are 
valid for Np ^ 0, otherwise there is no intraband contri- 
bution. Note that for high Fermi energy/low magnetic 
field, one recovers the familiar semi-classical Drude con- 
ductivity [Eqs. since the cyclotronic gap Aintra 
tends to the cyclotronic energy huc when many levels are 
occupied, TVf > 

Interband terms dominate at frequencies close or 
above the interband threshold, oj = 2Ep/h. These 
transitions involve the energy difference energy A„ = 
^/2hvp/lB{Vn+ 1 4- \/n), with n > Np. Its contribu- 
tion to the magneto-optical conductivity reads 



(inter) 



2e2 h^v^ 



h I 



B ^ — 



J2 (l + '5n,o)(2-(5„,Arp)x 



n=Nf 
1 



1 — ihuj/T 



(inter) 



2e2 h^vj. 



rA„ (1 - ihLj/T)^ + A2 /r2 

1 + 5np,o 



h 1%T^ il-thu/T)^ + A%jT^ 



(26) 
(27) 



with the notation. 



x{n) 



4e 



(29) 

(30) 
(31) 



The sign in the numerator in Eq. ([25]) must be chosen ac- 
cording to the requirements necessary to obtain a physi- 
cal solution (Appendix A). The function q{uj) then follows 
from Eq. (fT^ . 



(32) 



where er denotes the relative permittivity of the dielec- 
tric medium surrounding graphene. In order to proceed, 
we neglect the effect of interband transitions and approx- 
imate Eqs. (l24l) - ((25|) by its semi-classical analogue, [lOj 



e^2\Ep\ 



1 - ihn/r 



h r (i-ifti7/r)2 + A2/r2 



2Ep 



A/r 



h T (1 -iftf]/r)2 -I- A2/r2 



(33) 
(34) 



where A = heVpB /\Ep\ is the intraband cyclotron gap. 
The semi-classical expressions have the advantage of sim- 
plifying the notation and introducing the cyclotron en- 
ergy which is more used in the literature (albeit less ac- 
curate) than the intraband gap, Aintra- The crucial point 
to derive a compact expression for the dispersion rela- 
tion is to note that for typical frequencies ~TIIz and 
r ~ 0.01 eV, we have f x g ^ X"^, and therefore the 
expression for K(f2) can be approximated by 



Kin) 



x{n) 



(35) 



where we have chosen the appropriate sign in Eq. (1281) . 
(Note that the other solution has K{il) ~ and would 
correspond to a QTE mode.) Taking the square of both 
sides of Eq. ([55]). substituting the conductivity tensor 
components [Eqs. ([55 1) - ([M)) ] into X{Q), and employing a 
series expansion for small A/HD, and T/hfl, we obtain. 



q'" ^ er^ 



^^2 f2TTeMl\ 



\Epe^J ^ 



(36) 



VIII. APPENDIX C: SEMI-CLASSICAL 
DISPERSION 



We now derive an approximate formula for the semi- 
classical MPP dispersion. The first step is to solve Eq. ([2]) 
for K(r2); we obtain 



Kin) 



xin) ± ^4/(i7).9(f7) + A(r2)2 
mn) 



(28) 



where we have kept the terms up to second order in the 
small parameters, and assumed Ml/Ep ^ g^Z ~ 0.1 
(with go = 2e^ jh denoting the quantum of conductance). 
We remark that these approximations are consistent with 
the small wavelength limit, for which hn is typically 
larger than other energy scales. Assuming negligible 
damping g" ^ q', more precisely, requiring 



e2ft3f|3 



(37) 
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and taking q' ^ Q./c (non-retarded regime) in Eq. (15^ . 
we arrive at the final result, 

q'^^Vh'^'^2A^. (38) 

The magneto-plasmon spectrum [Eq. (|4])] follows imme- 
diately by expanding the latter expansion in the small 
parameter IS. /Ml. 
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